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A pioneering experiment [E. Schuster, E. Buks, M. Heiblum, D. Mahalu, V. Umansky, and Hadas 
Shtrikman, Nature (London) 385, 417 (1997)] reported the measurement of the transmission phase 
of an electron traversing a quantum dot and found the intriguing feature of a sudden phase drop 
in the conductance valleys. Based on the Friedel sum rule for a spinless effective one-dimensional 
system, it has been previously argued [H.-W. Lee, Phys. Rev. Lett. 82, 2358 (1999)] that the sudden 
phase drop should be accompanied by the vanishing of the transmission amplitude, or transmission 
zero. Here we address roles of strong electron-electron interactions on the electron transport through 
a two-level quantum dot where one level couples with the leads much more strongly than the other 
level does [P. G. Silvestrov and Y. Imry, Phys. Rev. Let. 85, 2565 (2000)]. We perform a 
perturbative conductance calculation with an explicit account of large charging energy and verify 
that the resulting conductance exhibits the transmission zero, in agreement with the analysis based 
on the Friedel sum rule. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

A pioneering experiment 0, on an Aharonov-Bohm 
ring containing a quantum dot reported the measure- 
ment of magnetic-flux-dependent interference signals and 
extracted the phase of the transmission amplitude of 
an electron traversing a quantum dot. The transmis- 
sion phase rises by tt as a new electron is introduced 
to the quantum dot and the measured phase profile 
near the Coulomb blockade resonance is in agreement 
with the Breit-Wigner formula The experiment re- 
vealed at the same time rather strange behaviors of the 
transmission phase; the transmission phase drops by tt 
almost suddenly in many conductance valleys and the 
transmission phase behaviors near two neighboring res- 
onance peaks are the same (up to 2tt) instead of show- 
ing a relative shift by tt. The sudden phase drop and 
in-phase resonances were later reproduced in other ex- 
periments 3. Recently it was reported @ that similar 
features persist in quantum dots with a relatively small 
number of electrons > 10. The experimental report 0,0 
induced a considerable amount of theoretical investiga- 
tions 0,1111 mis mill 

It is well known that the behavior of the transmission 
phase is constrained by the Friedel sum rule [l3] that re- 
lates the determinant of the scattering matrix with the 
number of electrons in the system. The constraint be- 
comes especially severe in one-dimensional systems and 
reduces to the form AQ/e = Aarg(<)/7r in strictly one- 
dimensional systems without any transverse degrees of 
freedom and without any side branches. Here Q is the 
total charge in a system, t is the transmission ampli- 
tude, and arg(t) represents the phase of t. Since AQ 
rises by e near each resonance peak but does not change 
in the conductance valleys, the strictly one-dimensional 
form of the Friedel sum rule predicts that the transmis- 
sion phase rises by tt only near each resonance. Thus in 



this prediction, a sudden phase drop in the conductance 
valleys is not possible and the transmission phase behav- 
iors near two neighboring resonance peaks should differ 
by TT. Therefore the experimental results in Ref. are 
not compatible with the strictly one-dimensional form of 
the Friedel sum rule. 

Shortly after the experimental report 0, three theo- 
retical calculations 0, iTol Hlf on the transport through 
a two-dimensional quantum dot or a quantum dot with 
a side branch were reported. It is found that when the 
quantum dot is linked to external electrodes via quantum 
point contacts so that a transport through the quantum 
dot becomes effectively one-dimensional, the transmis- 
sion phase does not follow the predictions of the strictly 
one-dimensional form of the Friedel sum rule and does 
show the sudden phase drops and in-phase resonances in 
close agreement with the experimental results 0. Inter- 
estingly the transmission amplitude is found to vanish 
(within the accuracy of the numerical calculation) when- 
ever the sudden phase drop occurs. The calculation re- 
sults are interpreted in terms of the Fano resonance |l4j , 
which can occur when the electron transport process me- 
diated by a level, which is weakly coupled to the elec- 
trodes, interferes with the electron transport process me- 
diated by a level, which is relatively strongly coupled 
to the electrodes. We remark that such fluctuations of 
the coupling strength require some deviations from strict 
one-dimensionality and are not possible in strictly one- 
dimensional systems. The Fano-resonance-based theory 
of the experimental results is further examined in the 
recent literatures |0, 0]. Unfortunately the electron- 
electron interaction effect is largely ignored in the Fano- 
resonance-based theory. Experimentally the appearance 
of the Fano resonance and its relation with the in-phase 
resonances were addressed in a quantum-dot-embedded 
Aharonov-Bohm interferometer |l7j . 

One of us examined the implications of the Friedel 
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sum rule in spinless effectively one-dimensional systems 
with time-reversal symmetry. It is found that the general 
Friedel sum rule AQ/e = [AlnDet(S)]/27r?, where S is 
the 2x2 scattering matrix ^T8!|, does not reduce to the 
strictly one-dimensional form since AlnDet(S)/2i can 
differ from Aarg(t) by tt whenever the transmission am- 
plitude vanishes identically (transmission zero). Similar 
conclusions are obtained in Ref. 19]. Thus the experi- 
mental results in Ref. do not violate the Friedel sum 
rule as long as the sudden phase drop is accompanied by 
transmission zero . Moreover it is argued based on 
the Friedel sum rule that the transmission zero is rather 
generic and robust in effectively one-dimensional systems 
since it is related to the topological structure ,21] of the 
Friedel sum rule. 

Silvestrov and Imry [T^ ] reported interesting effects 
of the strong electron-electron interaction on the trans- 
port through a quantum dot that is in the intermediate 
regime between integral and fully chaotic. They found 
that in such a semichaotic regime, certain single-particle 
levels in the quantum dot couple strongly with the leads 
while other levels couple very weakly with the leads [23] ■ 
Thus in such a semichaotic quantum dot, the widths of 
single-particle levels can differ from each other by orders 
of magnitude and in this sense, the situation addressed 
by Silvestrov and Imry is similar to that assumed for 
the Fano-resonance- based theory |i,[l3,[l3. But differ- 
ently from the Fano-resonance-based theory, they focused 
on the roles of the electron-electron interaction. In the 
strong interaction (or large charging energy) limit, they 
found that a new electron introduced to a quantum dot 
at a resonance peak may occupy a broad level with higher 
bare single-particle energy rather than narrow levels with 
lower bare single-particle energy, provided that the en- 
ergy gain via the hopping-induced downward level shift 
for the broad level can overcome the bare single-particle 
spacing. It was also found that as the gate voltage is 
increased towards the next Coulomb blockade peak, the 
energy gain disappears in the conductance valley between 
the two consecutive peaks and the electron in the broad 
level is transferred to an empty narrow level with lower 
bare single-particle energy. It was noted that the popu- 
lation switching of the broad level can be repeated over 
many consecutive resonance peaks, providing a natural 
explanation for the in-phase resonances. It was also ar- 
gued that the electron transfer from the broad level to 
a narrow level is responsible for the sudden phase drop. 
Thus the interaction-based theory of Silvestrov and Imry 
also produces a phase behavior in close resemblance with 
the experimental results 2]. 

In this paper, we revisit the population switching prob- 
lem addressed by Silvestrov and Imry We focus on 
the gate voltage region in the conductance valleys where 
the transmission phase drops suddenly by tt due to the 
electron transfer from the broad level to a narrow level. 
Our study is motivated by the observation that while 
the analysis ^ based on the Friedel sum rule predicts 
the sudden phase drop to be accompanied by a trans- 



mission zero and the prediction is indeed satisfied in the 
Fano-resonance- based theory H Elll [H , the relation 
between the sudden phase drop and the transmission 
zero is not clear in the interaction-based theory. Ref- 
erence Jl^ is rather focused on the equilibrium ground 
state configuration and does not address the transport 
properties in detail. A later publication l23 addressed 
the conductance behavior near the gate voltage range 
where the electron transfer occurs from the broad level 
to narrow levels. However it is rather focused on roles of 
the spin degrees of freedom and the relation between the 
sudden phase drop and the transmission zero is not ad- 
dressed. Recalling that the Friedel sum rule ^3 remains 
valid even in the presence of the electron-electron inter- 
action ^23] and that the transmission zero is an essential 
feature for the sudden phase drop to be compatible with 
the Friedel sum rule (at least when spin degrees of free- 
dom are not important), the relation between the sudden 
phase drop and the transmission zero in the interaction- 
based theory needs clarification. The goal of this paper 
is to verify this relation in the presence of the strong 
electron-electron interaction. 

The paper is organized as follows. The model Hamilto- 
nian and calculations of equilibrium and transport prop- 
erties are given in Sec. II. The results are discussed in 
Sec. III. 

II. MODEL AND CALCULATION 

For calculation, we use the following model Hamilto- 
nian, 

H = -ffdot + ^^int + ^^Icad + Ht, (1) 

where 

Hdot = eic[ci + snc^j^cn, (2) 

Hint = UcBclcic'j^CN , (3) 

i/iead = ^£(A:)4afc + Xl'^(0^J^/. (4) 

fe ; 

Ht =Y^ {tf-^clak + i^^cjvafc + iJ.c.) (5) 
fe 

+ J2{t^i'^<^\b,+tfc^^b,+H.c.), 
I 

where , 5/ , and Ci {i — I, N) denote the electron an- 
nihilation operators for electrons in the left and right 
electrodes and in the quantum dot, respectively. Uqb 
is the constant charging energy due to electron inter- 
action in a quantum dot. tp-* is a hopping coefficient 
of electrons moving from the lead j (L for the left lead 
and R for the right lead) to the dot state i. In order to 
make the model as simple as possible while keeping the 
main physics of the broad level population switching in 
a semichaotic dot, the quantum dot is assumed to have 
only two single-particle levels (levels i = 1 and i = N), 
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where one of them (level N) is a broad level with larger 
couplings with the leads and the other (level 1) is a nar- 
row level with smaller couplings with the leads, ei and 
En denote single-particle energies of dot states |1) and 
|A''), respectively, which can be modulated by the gate 
voltage Vg, 



(6) 



where — e(< 0) is the electron charge and k is a dimen- 
sionless constant depending on the gate geometry. e(k), 
the energy of an electron in the leads, is defined as 



e{k) 
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2m 
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where ep is the Fermi energy. Note that the one- 
dimensional dispersion relation is assumed for the leads 
in order to describe a quantum dot coupled to external 
electrodes via quantum point contacts. We are interested 
in temperature regimes well above the Kondo tempera- 
ture and for simplicity, electron spin degrees of freedom 
are ignored. This model is the exactly same as the model 
Hamiltonian analyzed in Ref. |12j except for a trivial gen- 
eralization from a single lead to two leads. Below electron 
transport from one lead to the other lead via the quan- 
tum dot will be investigated. 

Reference addressed the situation (though only 
one lead is taken into account) where the broad level lies 
above the narrow level {en > ^i)- By using the nonde- 
generate perturbation theory with H^ot + Hmt + i?icad 
as an unperturbed part and Ht as a perturbation, it 
was noted that the large second order energy correction 
(oc I^TvP) for the broad level can be larger than the bare 
energy level splitting ejv — ei in the semichaotic dot. Af- 
ter careful comparison of the two possible configuration 
of the dot state in the gate voltage range where only one 
electron is allowed in the dot, it was concluded that it is 
energetically favorable for an electron in the dot to oc- 
cupy the broad level in the lower half of the gate voltage 
range and to occupy the narrow level in the upper half 
of the gate voltage range. In the next two subsections, 
we first reexamine the problem of the ground configura- 
tion in the gate voltage range where the electron transfer 
occurs from the broad level to the narrow level. 



a single electron in the dot. For a given energy eigenket 
of H with = E\^), its projection Pil^") onto 

the subspace with only one electron in the dot becomes 
an energy eigenket of i?eff , 



Hii + Hi 



1 



E-Hn 



-Hqi + Hi2 



1 



E — H22 



H21, (8) 



where H,nn = PmHPn and P„ is a projection operator 
onto the subspace of states with n electrons in the dot. 
Equation (jSJ can be derived by using the Schrieffer- Wolff 
transformation 25]. Note that the resulting Schrodinger 
equation i?off [^'il^')] = -^[-Pil^)] is a self-consistency 
equation in the sense that Hcs itself contains the exact 
energy eigenvalue E. As long as the exact E is used, the 
transformation from H to HcS is exact. 

Next we apply the perturbation theory to Hcs- We de- 
compose Hes into the unperturbed Hamiltonian iJunpert 
and the perturbation H' as follows: 



H, 



cff 



(9) 



unpcrt 



H 



Pi Eavc (c\ci + cjyCAr^ Pi 
+Pl [Hint + Hiead] Pi , 

ei - En f ^ t 

Pi 2 V"^!*^! ~ ^N^^N 



+H10 



E- Hi 



-Hoi + H12 



Pi 
1 



00 



E — H22 



-H 



21, 



where Save = (£i+£Ar)/2. Note that H^^t is split into two 
parts, the average energy part (the first term in i/unport 
proportional to Eave) and the level spacing part (the first 
term in H' proportional to £jv — ei). Since the level spac- 
ing part is included in H' together with the last two terms 
of H', which are of order and responsible for the en- 
ergy correction by hopping, the perturbation calculation 
with H' as a perturbation provides a desired equal foot- 
ing treatment of the two competing energy scales. For 
the purpose of the first order perturbation in H' , E in 
H' may be replaced by i?'"-* (unperturbed energy eigen- 
value of iJunport) since the difference E — E'^^^ affects only 
higher order perturbation calculations. 



Ground state 



A. Effective Hamiltonian 

In order to systematically address the gate voltage 
range very close to the electron transfer, where the bare 
energy level spacing en —Ei competes with the second or- 
der energy correction via hopping, an equal footing treat- 
ment of the two competing energy scales is desired. For 
this purpose, we derive in this subsection an effective 
Hamiltonian that facilitates such an equal footing treat- 
ment. Recalling that the electron number fluctuation in 
the dot is strongly suppressed in the conductance valley 
under investigation, one can construct an effective Hamil- 
tonian that acts only on a special subspace of states with 



Unperturbed ground states of -ffunpcrt are doubly de- 
generate and given by 

\ipi) — cjjvacuum), 
Wn) = c]y I vacuum), 

where | vacuum) is defined as the state with zero electron 
in the dot and single-particle-levels in the leads com- 
pletely filled up to the Fermi level. The correspond- 
ing unperturbed ground state eigenenergy is i?Q°^ = 
Eavc + Efc<fep £{k) + Y.i<kF £(0- The degeneracy is hfted 
by the perturbation H' . According to the degenerate per- 
turbation theory, the first-order energy corrections are 
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eigenvalues of H' within the subspace spanned by 
and \tfM)'- 



\H')ni {H')nn 



{ipN\H'[E^E, 



AO) 

^0 
p(0) 



\ipi){^N\H' (S: 



E, 



(0) 



E, 



(0) 



I \V^n) 



(10) 



where 
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Here T 



(L) 



27r 



.(R) 



C/CB 

£(fc)=Ol 

(L) 



and A6l(L) = 
and are de- 



fined in a similar way. The correct zeroth-order ground 

state I ■00°^) = + P\^n) is also determined by 

Eq. Hl()|l since (^) is a normalized eigenvector that cor- 
responds to the lower eigenvalue of Eq. . Note that 
when the off-diagonal matrix elements of Eq. 1)10(1 are 
ignored, the eigenvectors are (g) and and the first 
order energy corrections, which are the diagonal matrix 
elements of Eq. (fTUIl . reproduce the result in Ref. 
except for trivial corrections due to the extension from a 
single lead to two leads. The dashed lines in Fig. Qla) 
depict the energy eigenvalues of Eq. pO|) without the off- 
diagonal matrix elements. Note that two levels cross at a 
certain gate voltage, at which the electron configuration 
changes suddenly from ("j*) to (J). The inclusion of the 
off-diagonal matrix elements modifies the result near the 
level crossing point. As Fig. Qla) shows, the level cross- 
ing is avoided (solid line), and the configuration change 
of the dot level occurs smoothly over a finite gate voltage 
range [Fig. ^b)]. Later it turns out that the deviation 
of the eigenvector (^) from (J) and (°) near the avoided 
crossing is important for the relationship between the 
sudden phase drop and the transmission zero. A similar 
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FIG. 1: (Color online) (a) Energy eigenvalues of Eq. IIUII as 
a function of eavo/f/cB. When the off-diagonal matrix ele- 
ments are ignored, two energy levels cross (blue short-dashed 
line and red long-dashed line). When ofF-diagonal matrix el- 
ements are taken into account, the energy level crossing is 
avoided (solid lines), (b) (blue short-dashed line) and 
|/3p (red long-dashed line) as a function of eavo/C/cB, where 
(^) and (jf^Q.) are two eigenvectors of Eq. IIOII corresponding 
to lower and higher energy eigenvalues shown in (a). Con- 
ductance (black solid line) is also shown. The conductance is 
multiplied by a proper factor to fit in the graph box. Thus 
the marked conductance values are in arbitrary units. Note 
that the conductance shows a dip near the avoided crossing. 
In this plot, the time reversal symmetry is assumed and the 
conductance becomes zero (transmission zero) at the center 



of the dip. Here t^^V^f = b CV^i"' = 10, t^Vt^"' = 11. 
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avoided crossing is reported in Ref. l2a| which addresses 
roles of the spin degrees of freedom |2^ . 



C. Conductance 

A standard linear response calculation in the Appendix 
A results in the following formula for the linear response 
conductance. 



G=— hm {M{E, 



Huj)huj] (^'o|(57Vdiff|5'„) 



(11) 



where Af{E) is the many-body density of states, |^„} is 
the n-th excited states of H with energy En, SlSf^^^ = 
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FIG. 2: Schematic plots of four groups of excited states Iv&n) 
that contain an electron-hole pair excitation in the leads with 
the excitation energy hui. 



Next we evaluate («'o|<5iV'^''^|*„). Let \<po) = Pi\^o), 
\(pn) = Pil^'n), which are solutions of HcS (Eq) \ipo) = 
Eq\ipq) and Hcs (En) \ipn) = En\ipn)- According to the 
SchriefFer- Wolff transformation, \'^o) and are re- 

lated to \ipo) and !</?„) as follows, 

1^0)= ^ ^ rr Hoi\ipQ) + \ipo) + —^-—H2l\'fo), 



Eq — H{ 



00 



Eq — H22 



(12) 



|^«)= ^ HoiWn) + \Vn) + „ H2l\(pn) 



En — Hn, 



En — Ho 



(13) 



where the first, second, third terms represent the projec- 
tions of I'I'o) and \'^n) onto the subspace with zero, one, 
two electrons in the dot, respectively. By using Eqs. I|12|) 
and lO and recalling that PsSN'^'^Pt = for s ^ t 
(s, t = 0, 1, 2), one finds 



represents the average of 



and (^'o|(5Ardiff]^'„} 
2 

(^'o|(5iV'*'**|«'„) over states l^-,,) with E,, = 

Eo + huj- Note that Eq. Hll(l is expressed in terms of the 
energy eigenstates and eigenvalues of the original Hamil- 
tonian H . 

the limit uj 0+, M(Er) + Hlu) and 
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^ 0+, J\f{Eo + hio) 

for the original Hamiltonian 

H can be evaluated by using the effective Hamiltonian 
Hes- Here we evaluate Eq. ((TT|l up to the second or- 
der in r and En — £i- We first evaluate M {Eq + hiu), 
which is governed, in the w — > 0^ limit, by low energy 
excitations. Low energy excitations of the system are 
electron-hole pair excitations in the leads and excitations 
in the dot configuration. But as demonstrated in Sec. II 
B, the excitation in the dot configuration has a finite 
excitation energy due to the avoided crossing and thus 
excitations in the dot can be neglected in the w ^ 0+ 
limit On the other hand, excitations in the leads 

can have infinitesimal excitation energies and thus con- 
tribute to JV {Eq + tujj) . Figure Q shows schematically 
four groups of excited states with a single electron- 
hole pair excitation in the leads. States IvP^) in groups 
(a) and (b) may be neglected for the density of states 
evaluation since |(5'o|(SiV'^''^|^'„)p = up to second or- 
der in r and En — £i- On the other hand, for states 
in groups (c) and (d), K^ol^iV^iff |*„)|2 _ q(y^2 

and groups (c) and (d) contribute hujAfsp\eF)J^sp\eF) 
to Af{Eo + hu), where Af^p\eF) and Af^p\eF) are sin- 
gle particle densities of states of the left and right leads, 
respectively, at the Fermi energy. We remark that states 
l^'n) with n pairs of electron-hole excitations also con- 
tribute to Af {Eq + Hlu) but their contribution scales as 
cj^n-i^ which is negligible in the w — > 0"*" limit. 



(vI/o|^7V^'"|vI'„) 



((po\Hw 
+ (■^01' 

+ (',^'01-^12 
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Eq — Hoo 
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6N' 



diff 
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En — H; 



■SN- 



En — Hqo 
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diff 



22 



En — H; 



Hoi\(pn) 



-H2l\(pn 



(14) 



22 



Note that {'i> o\S N'^^^ „) is now expressed in terms of 
\ipo) and \(pn), which are eigenkets of -ffeff. We evalu- 
ate the three terms on the right hand side of Eq. H14|) 
up to the first order in F and En — Ei by performing 
the perturbation theory calculation for Heg. After some 
calculations, it can be verified that the first and third 
terms on the right hand side of Eq. H14|l are finite in 
the w ^ O"*" limit while the second term is proportional 
to uj-\ Thus in the w ^ 0+ limit, (^'o|(5iV'i'^|*„} 
completely governed by the projections of 1^*0 ) and 
onto the subspace with one electron in the dot. Note that 



the 



dependence of 



is canceled by 



[J^ {Eq + hjj) hjj] (X w^, producing a finite G [Eq. l(TT1l ] 
in the cj ^ 0+ limit. Combined with the result for 
TV [Ea + hid), one obtains 



G 



27re" 



1 



Eave + C^CB 



-A_B* 



(15) 



where 
A 



+ 



+ (L)* 



a*t 



S_ = PH^^' - a*t' 



(R)' 



Here (^) and {_^,) are the normalized eigenvectors of 
Eq. H1U|I corresponding to its lower and higher eigenval- 
ues, respectively. 
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Before we discuss in the next subsection the imphca- 
tions of Eq. (|15|l regarding the transmission zero, we dis- 
cuss the physical interpretation of Eq. (|15|l in terms of 
elastic cotunnehng processes. For this purpose, we first 
introduce two linear combinations |D1) and |D2) of the 
dot states |1) and |A^), 



(a) 



(b) 



|D1) =a|l)+/3|iV), 
|D2) = /3*|1) -a*\N). 



(16) 
(17) 



Since (^) and (_^.) are the two normahzed eigenvec- 
tors of Eq. (|HJII . |D1} and |D2) represent two effective 
single-particle eigenstatcs in the dot. In terms of |D1) 
and |D2), Eq. IjlSI) has a simple physical meaning; A+ 
and B+ represent effective hopping matrix elements from 
|D1) to the left and right electrodes, respectively, and A- 
and B- represent effective hopping matrix elements from 
|D2) to the left and right electrodes, respectively. In this 
description, the first term — ■^+B'!^ in Eq. (|15() can be 
interpreted as the amplitude of the cotunnehng process 
in Fig. ini^a); an electron in |D1) hops to the left elec- 
trode first (v4+) and then another electron in the right 
electrode hops to the |D1) (-B+). The energy of the in- 
termediate virtual state with respect to the initial state 
is — £ave(> 0). Similarly the second term - — \u^^ A^B*_ 
in Eq. H15|l can be interpreted as the amplitude of the 
cotunnehng process in Fig. I^b) ; an electron in the right 
electrode hops to the empty level |D2) {B*J) and then the 
electron in |D2) hops to the left electrode {A_). The en- 
ergy of the intermediate virtual state with respect to the 
initial state is Eavo + C^cb(> 0). Note that both cotunnel- 
ing processes are elastic since the initial and final states 
have the same energy. To obtain the total transition 
amplitude, the transition amplitudes for the two elastic 
cotunnehng processes should be summed up coherently. 
A careful comparison of the two final states in Fig. ISJa) 
and (b) indicates that they differ by a pair-wise electron 
exchange of the electron in |D1) and an electron in the 
left lead. This pair- wise exchange gives rise to a relative 
phase factor (—1), which explains the (— ) sign in front 
of ^ — ^+-8+ in Eq. H15|l . In this description in terms 
of the effective dot states |D1) and |D2), Eq. H15|l is also 
consistent with the (second order) golden rule formula. 
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p{Ef) 



E f—Ei 



(18) 



where the effective perturbation H represents the hop- 
ping between the leads and the effective dot states |D1) 
and |D2), and w^i^^y] is related to G via ew\^i\^^f^ ~ 

GAF a.xidp{Ef) \e,^e. = e^VN^l\eF)M^f\eF)- Here 
Ay denotes the voltage difference between the two leads. 



D. 



Transmission zero 



We examine in this subsection the Eavc-dependence of 
G and demonstrate the appearance of the transmission 
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FIG. 3: Two elastic cotunneling processes of two electrons, 
(a) and (b) describe hopping via the effective dot states jDl) 
and |D2), respectively. The arrows indicate the electron hop- 
ping direction. The number (1) and (2) above the arrows is to 
denote which hopping occurs first (1) or next (2). Associated 
hopping amplitudes are also denoted. 



zero. Experimentally, Save can be modulated by the gate 
voltage applied to the dot [Eq. ©]. To get an insight, we 
first examine a special case of symmetric hopping with 
t^i^ = 4^"* = ti and = — t^- In this case, both 
A+B^ and A^BZ_ are real and thus the total amplitude 
(within I • • ■ I) in Eq. (|15() is real. Based on behaviors of 
a and /3 in terms of s^vc, which was reviewed in Sec. 11 
B, one can estimate the conductance in the following two 



regions, (1) -1 < ^ < -1 and (2) -i < < 0. In 
case (1), a ~ 1 and /3 ~ 0, and the total amplitude of 
the two cotunneling processes in Eq. (|15|l becomes 



\ti\' 



+ 



\tN\' 



(19) 



which is a positive real number as shown in Fig. ^ In 
case (2), a ~ and /3 ~ 1, and the total amplitude in 
the region (2) reduces similarly to 



\tN\' 



\h\' 



Ucb' 



(20) 



fr^ = — i, and the 

Ucb 2 ' 

Note that the transmission 



which is a negative real number as shown in Fig. 0] Since 
the total amplitude varies continuously and remains real 
in the whole range — 1 < < 0, Eqs. (fll^ and in- 
dicate that the total amplitude should vanish identically 
at a particular value of £avc near 
transmission zero appears 
zero appears in the region where the ground state dot 
configuration changes. 

Next we examine a generic situation with non- 
symmetric hopping and demonstrate that the transmis- 
sion zero persists even in this situation if the system has 
a time reversal symmetry. In a time-reversal symmetric 
system, all hopping matrix elements t^'"^ , t^^"* , t j^"* , t^"* 
can be taken real upon proper gauge transformations of 
the electron annihilation operators, a e^^^a, b — > e*^''5, 
c — > e'^^c. Then Eq. UlUI) becomes a real symmetric ma- 
trix, and both a and (3 and also A+B'^, A-B*_ become 
real. Hence the total amplitude [expression within the 
absolute value symbol in Eq. H15|) ] again becomes real 
and it can be verified that the amplitude changes its sign 



FIG. 4: Schematic plot of the transmission amphtude in two 
regions: (1) -1 < ^ < _i and (2) -i < |^ < 0. The 
regime around ~ ~ 5 requires a numerical evaluation [see 
Figs.IHb) and El'?'' 



similarly to Fig. 0] signaling the occurrence of the trans- 
mission zero. Figure Q^b) shows numerical calculation 
of conductance for a time-reversal symmetric system in 
the dot configuration crossover regime Eavc/C^CB ~ ~l/2. 
Note that the transmission zero happens indeed when 
the configuration change occurs. We have assigned cou- 
pling coefficients as follows: tj'^-'/^^'' = l,t7v''/4^'' ~ 

10,i^V4^^ = 11- 

Next we consider systems without the time-reversal 
symmetry. Then not all hopping matrix elements can 

(7) 

be taken real, and the phases of t] do affect G. From 
Eqs. lfTn|l and l(T5|) . it can be verified that the phase 
dependence of G occurs only via the total phase of 

^'"^^TV^ ^jv^"'^!^'' instead of individual phases. Figure 
shows the G vs. £avc/t^CB graphs as a function of the to- 
tal phase. Note that the transmission zero appears only 

when t1'\^^'^ ^Tv^''^!^^ is real (total phase of 0, tt, or 27r), 
and the transmission zero is replaced by a small but fi- 
nite conductance dip when t'^^^t^^'^ t^N ^^t'f^^ is not real. 

When t^^^tj^^%^^4^^* is real, the all tp'^'s can be taken 
real via proper gauge transformations of the electron an- 
nihilation operators, implying that the system has the 
time-reversal symmetry. Thus this total phase depen- 
dence of the transmission zero illustrates the role of the 
time-reversal symmetry for the transmission zero. This 
result is consistent with the analysis based on the 
Friedel sum rule and also agrees with the result for 
the Fano-resonance-based theory. 

E. Transmission phase 

Here we consider the phase of the transmission ampli- 
tude through the quantum dot. Experimentally testable 
definition of the transmission phase requires an interfer- 
ence configuration such as the Aharonov-Bohm ring used 
m Ref. 0. For this purpose, we introduce a direct hop- 
ping channel from the left to right leads, not mediated 



-0.4934 -0.4932 -0.493 -0.4928 -0.4926 

^ave /UcB 



(b) Total phase: tt to 2/1 




-0.4934 -0.4932 -0.493 -0.4928 -0.4926 

'^ave /UcB 

FIG. 5: (Color online) Evolution of the conductance 
VS . £ ave /UcB curve as a function of the total phase of 

^I'^'^jv' ^iv^'^i'^^ • (8-) From the total phase of (red, labeled 
as 0) to TT (blue, labeled as tt) in steps of 7r/8 and (b) from n 
(blue, labeled as tt) to 2n (red, labeled as 2tt) in steps of 7r/8. 
The transmission zero occurs only when the total phase is 0, 
TT, or 2n. Note that the horizontal axis is magnified consid- 
erably compared to that of Fig. ^b), so that the whole span 
of the horizontal axis in these figures covers a narrow portion 
of the Eavo/f/cB range near the conductance dip in Fig. 0^b). 
The units for the conductance in (a) and (b) are the same as 
that used in Fig. 0^b) . 



by the dot levels. Equation (|15fl is then modified to 

1 1 ^ 

X A^l + A_B*_+t^^^^^^'> , 

S^avc Save + UcB - 

(21) 

where is the transmission amplitude of the 
direct hopping from the left to right leads. Thus 
the phase of the transmission amplitude through a 
quantum dot is nothing but the phase of j-^A+B* -|- 
:pTj — A^B* . Figure El shows the transmission phase 
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(a) Total phase: to tt 
5 7T F I ^= 




-0.495 -0.494 -0.493 -0.492 -0.491 

£ave/UcB 



(b) Total phase: tt to 2/1 
5 TT h =i 




-0.495 -0.494 -0.493 -0.492 -0.491 

Save/UcB 

FIG. 6: (Color online) The transmission phase vs eave/(7cB as 
a function of the total phase of t^^H^^^ i^'t''^' . (a) From the 
total phase (red line, labeled as 0) to tt (blue line, labeled 
as tt) and (b) from tt (red line, labeled as tt) to 2n (blue line, 
labeled as 27r). For clarity, the curves are shifted vertically 
by proper amounts. 



versus Eave/t^CB graphs as a function of the total phase 

of tf^H'"^^ t^^^t'f^^ . In systems with time-reversal sym- 
metry (with the total phase 0, tt, or 27r), the transmission 
phase changes abruptly by tt at transmission zero even 
though the dot configuration of the ground state changes 
smoothly. In systems without time-reversal symmetry, 
the transmission phase changes smoothly by tt. Note that 
when the total phase is between tt and 2tt, the transmis- 
sion phase drops by tt as £avc/C^CB is increased, and when 
the total phase is between and tt, it rises by tt. 



III. CONCLUSION AND DISCUSSION 

In the interaction-based theory [T^ of the experimen- 
tal results in Ref. the dot configuration change in 
conductance valleys is responsible for the sudden phase 



drop. We performed the perturbative conductance cal- 
culation up to the first nonvanishing order in the elec- 
tron hopping matrix elements (fourth order in tp'') and 
found that the dot configuration change is accompanied 
by the transmission zero. According to the Friedel sum 
rule, transmission zero is a mandatory feature for the 
sudden phase drop in a spinless time-reversal symmetric 
one-dimensional systems ip], and our calculation shows 
that the interaction-based theory 0| is consistent with 
the prediction of the Friedel sum rule (at least up to the 

fourth order in ip''). 

A few remarks are in order. We first compare the 
interaction-based theory and the Fano-resonance- 
based theory [3>E)I3 of the sudden phase drop and the 
in-phase resonance. The two theories are similar in the 
sense that they both assume the coexistence of strongly 
coupled levels and weakly coupled levels. But regarding 
the roles of electron-electron interaction, the two theories 
are different; While the interaction is essentially ignored 
in the Fano-resonance-based theory 9, 10, llj, it plays 
a crucial role in the interaction-based theory One 
manifestation of this difference is the relevance or irrel- 
evance of the signs of the hopping matrix elements. For 
simplicity of illustration, we assume a time-reversal sym- 
metric system, where all hopping matrix elements may 
be assumed to be real. According to Ref. [l^, which 
addresses the role of the hopping matrix signs in the 
context of the Fano-resonance-based theory, the sudden 
phase drop (and also transmission zero) may or may not 
appear between two consecutive Coulomb blockade reso- 
nance peaks, depending on the relative signs of the hop- 
ping matrix elements. In the interaction-based theory, 
in contrast, the signs of the hopping matrix elements are 
largely irrelevant to the existence of the sudden phase 
drop (and also transmission zero). As demonstrated in 
Figs. and ini the sudden phase drop (and transmission 

zero) persists even when the sign of one of tp'''s changes 

(and thus the total phase of tf^^'t'"^'' t^^hf^'^ is altered 
from to TT or from tt to 0). Only its position is shifted by 
a small amount (fcavc/C^CB ^ 1) due to the sign change. 

Our discussion so far neglected the spin degrees of free- 
dom. When the transport is spin-dependent or when 
spin-flip scattering becomes possible, it is crucial to take 
into account the spin degrees of freedom. In this case, 
the electron transmission from one electrode to the other 
is described by four (instead of one) transmission am- 
plitudes, i|f , til, ^Ti' ^iT' and each transmission ampli- 
tude may have different gate voltage dependence. Refer- 
ence |23| addressed spin effects of the interaction-based 
theory. In the limit t'"^^^^ ^ t'"^^^\ the total conduc- 
tance (cx |i|tP -I- l^uP + I^TiP + kiTp) shows interest- 
ing features (see Fig. 3 in Ref. [23) near the gate volt- 
age for the dot configuration change. Unfortunately the 
individual behaviors of the four amplitudes are not ad- 
dressed. Further insights into the experimental results 
in Ref. 2| may be obtained from spin resolved studies 
of the transmission amplitudes and their dependence on 
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time-reversal symmetry. 
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APPENDIX A: LINEAR RESPONSE 

To obtain the hnear response conductance of the sys- 
tem described by the Hamiltonian H , we first introduce 
the voltage bias into the total Hamiltonian, 



k 



(Al) 



and use the method of the adiabatic turning on, 

H-tot{t) = H + -F(i)i?bias, 



where F{t) — e^* cos cut with 77 ^ 0+ and the limit 
a; ^ 0"^ will be taken to obtain the DC conductance [s^l ■ 
From the Kubo formula [sij . the current expectation 
value I{t) in the linear response regime is readily ob- 
tained as follows: 



lit) 



1 /■* 

hm — / dt'F{t ) 



bias 



(A2) 



where |^'o) is the exact ground state of H and / is the 
current operator given by 



dt \ 2 



2ih 



.(A3) 



Here iV^^) ee Y.k<^W [N''^'^ = Y.k^W) counts the 
number of electrons in the left (right) lead, and Q*-'"-' = 
-e7V('"), Qf^) = -e7V(^). After inserting the closure 
X]n=o i^^^o abovc cxprcssiou, where |5'„) is 

the n-th excited states of H with energy £"„, one per- 
forms the integration over time t and take 77 — > 0"*" limit 
to obtain the expression for the conductance G, 



G- 



m 

= 1™ ^V<5(i?„-i?o-M^^I(*ol^'''''l^«)l', 



(A4) 

where N'^'^ = A^^^) - iV^^l We use the identity 
(^fo|iVdiff|,Ifg) and obtain the final expression, 



G^— lim [M (Eo + huj) huj] (*o|5iVdiff|4r \ 



(A5) 



where M [E) is the many-body density of states 
and 



(*o|(57Vdiff'|^J 
2 



represents the average of 
over states with En — Eo + hw. 
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